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Abstract 

We present the lattice gauge theory approach to evaluating non-perturbative hadronic 
interactions from first principles. We discuss applications to glueballs, inter-quark poten- 
tials, the running coupling constant, the light hadron spectrum and the pseudoscalar decay 
constant fp. 



1 Introduction 

The theory of the strong interactions is accurately provided by Quantum Chromodynamics 
(QCD). The theory is defined in terms of elementary components: quarks and gluons. The only 
free parameters are the quark mass values (apart from an overall energy scale imposed by the 
need to regulate the theory). This formulation is essentially the unique candidate for the theory 
of the strong interactions. The only feasible way to describe a departure from QCD would be 
in terms of quark and gluon substructure. At least at the energy scales up to which it is tested, 
there is no evidence for such substructure. 

QCD provides a big challenge to theoretical physicists. It is defined in terms of quarks and 
gluons but the physical particles are composites: the mesons and baryons. Any complete descrip- 
tion must then yield these bound states: this requires a non-perturbative approach. One can see 
the limitations of a perturbative approach by considering the vacuum: this will be approximated 
in perturbation theory as basically empty with rare quark or gluon loop fluctuations. Such a 
description will allow quarks and gluons to propagate essentially freely which is not the case 
experimentally. The true (non-perturbative) vacuum can be better thought of as a disordered 
medium with whirlpools of colour on different scales. Such a non-perturbative treatment then 
has the possibility to explain why quarks and gluons do not propagate (ie quark confinement). 
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The main contender for a non-perturbative description of QCD is the lattice gauge theory 
approach. Following the ideas of Wilson, space-time is replaced by a discrete grid (the lattice) 
but gauge invariance is retained exactly. Using also periodic boundary conditions in space and 
time, the system will have a finite number of degrees of freedom: the gluon and quark fields 
at the lattice sites. Actually gauge invariance implies that the gluon field lies on the links of 
the lattice (they correspond to the mechanism that allows colour orientation at different sites to 
be compared). This finite number of degrees of freedom implies that the theory is a quantum 
many-body problem rather than a field theory. 

The step which makes this many-body problem tractable is to consider Euclidean time. 
This is perhaps the step which is the most difficult to assess. Formally the quantities to be 
evaluated can be expressed as Green functions: vacuum expectation values of products of fields. 
The existence and properties of these Green functions under a Wick rotation to Euclidean time 
are widely used in perturbative treatments. A non-perturbative determination of the Green 
functions in the Euclidean time case will lead directly to quantities of relevance in the physical 
Minkowski time case (eg masses and matrix elements). 

With the Euclidean time approach, the formulation of QCD (consider, for example, the 
functional integral over the gauge fields) is converted into a multiple integral which is well 

defined mathematically. For a lattice of sites with a colour gauge group of SU(3), this would 
be a 8 X 4 X dimensional integral (8 from the colour group manifold, 4L'^ from the gauge fields 
on each link). For any reasonable value of L, this is a very high dimension indeed. Simpson's 
rule is not the way forward! The standard approach is to use a Monte Carlo approximation to 
the integrand. This is implemented in an "importance sampling" version so that a stochastic 
estimate of the integral is made from a finite number of samples (called configurations) of equal 
weight. The construction of efficient algorithms to achieve this is a topic in itself, ffere 1 
will concentrate on the analysis of the outcome, assuming that such configurations have been 
generated. 

So what we have at our disposal is a set of samples of the vacuum. It is then straightforward 

to evaluate the average of various products of fields over these samples - this gives the Green 
functions by definition. The Green function can then be continued from Euclidean to Minkowski 
time (in most cases this is trivial) and compared to experiment. 

The validation of the lattice approach calls for a series of checks that everything is under 
control. 

• the lattice spacing should be small enough (discretisation errors) 

• the lattice must be big enough in space and time (finite size errors) 

• the statistical errors must be under control 

• Green functions must be extracted with no contamination (eg a ground state mass could 
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be contaminated with a piece coming from an excited state) 

• both the quark contribution to the vacuum (sea quarks) and the quark constituents of 
hadrons (valence quarks) are usuahy treated by using larger mass values than the experi- 
mental ones and then extrapolating. This extrapolation must be treated accurately. 

• a particularly severe approximation is to treat the sea-quarks as of infinite mass. This 
corresponds to neglecting quark loops in the vacuum and is computationally very much 
faster. This is known as the "quenched approximation". 



The most subtle of these is the discretisation error. In order to extract the continuum limit 
of the lattice, one must show that the physical results will not change if the lattice spacing is 
decreased further. This is subtle because the lattice spacing is not known directly - in effect, 
it is measured. The lattice simulation is undertaken at a value of a parameter conventionally 
called p. In the limit of small coupling g'^, where perturbation theory applies, /3 = Thus 
large f3 corresponds to small g^. Now, perturbatively, the coupling g^ corresponds to the lattice 
spacing a as 

g f« — 2 where bo = — — 

Oq log a ^ Ibvr^ 

for N ^ flavours of quarks. For the case of interest, Nj < 3, this corresponds to small values of 
g^ at small distance scale a - as expected from asymptotic freedom. 

The perturbative argument is appropriate to the study of results at large /? (small g'^). 
We will find that lattice simulation of QCD uses values of /3 ~ 6 and hence bare couplings 
corresponding to a = g"^ /{An) ~ 0.08. In the pioneering years of lattice work, this was thought 
to be a sufficiently small number that the perturbation series would converge rapidly. One of 
the major advances, in recent years, has been the realisation that the bare lattice coupling (our 
g^ above) is a very poor expansion parameter and the perturbation series in the bare coupling 
does not converge well at the /? values of interest. The theoretical explanation for this poor 
convergence is that the lattice Lagrangian differs from the continuum Lagrangian and allows 
extra interactions. These include tadpole diagrams which have the property that they sum 
up to give a contribution that involves high order terms in the perturbation series. The way to 
avoid this problem with tadpole terms is to use a perturbation series in terms of a renormalised 
coupling - rather than the bare lattice coupling. I return to this topic when discussing the lattice 
determination of as later. 

This change of attitude to the method of determining a from /? has had considerable implica- 
tions for lattice predictions, as I now explain, since we wish to work in a region of lattice spacing 
a where perturbation theory in the bare coupling is not precise. Because of this, in practice, a 
is determined from the non-perturbative lattice results themselves. Thus if the energy of some 
particle is measured on the lattice, it will be available as the dimensionless combination E. 
Prom a value for E in physical units, then a can be determined since, on dimensional grounds, 
E = Ea. Furthermore, by increasing /3, the change in the observable E gives information about 
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the change in a since E is fixed - assuming it is the physical mass. This should allow a calibration 
of P in terms of a to be established. 

This procedure is overly optimistic, however. The discretised lattice theory is different from 

the continuum theory on scales of the order of the lattice spacing a. For the Wilson action 
formulation of gauge theory, this implies that the continuum energy Ec is related to the lattice 
observable E{P) = Ea as 

^ = E, + 0{a') 

A direct consequence is that the ratio of two energies (of different particles, for example) will 
have discretisation errors of order a?. 

Thus, to cope with discretisation errors, the procedure required is to evaluate dimensionless 
ratios of quantities of physical interest at a range of values of the lattice spacing a and then 
extrapolate the ratio to the continuum limit (a — >^ 0) . 

Note that when the fermionic terms are included, the discretisation error is of order a for 
the Wilson fermionic action. By adding further terms in the fermionic action, the error can be 
reduced - to order aa for the SW-clover fermion formulation. 



2 Glueball Masses 



I choose to illustrate the workings of the lattice method by describing the determination of 
the glueball spectrum. Of course, glueballs are only defined unambiguously in the quenched 
approximation - where quark loops in the vacuum are ignored. In this approximation, glueballs 
are stable and do not mix with quark - antiquark mesons. This approximation is very easy to 
implement in lattice studies: the full gluonic action is used but no quark terms are included. 
This corresponds to a full non-perturbative treatment of the gluonic degrees of freedom in the 
vacuum. Such a treatment goes much further than models such as the bag model. 

The glueball mass can be measured on a lattice through evaluating the correlation C{t) of 
two closed colour loops (called Wilson loops) at separation t lattice spacings. Formally 

C{t) =< 0|G(0)G't(t)|0 >=J2^^< 9^m9^it) > < = ^ hl^^-"^'' 

1=0 i=0 

where G represents the closed colour loop which can be thought of as creating a glueball state 
gi from the vacuum. Summing over a complete set of such glueball states (strictly these are 

eigenstates of the lattice transfer matrix where exp —rhi is the lattice eigenvalue corresponding 
to a step of one lattice spacing in time) then yields the above expression. As t ^ oo, the lightest 
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glueball mass will dominate. This can be expressed as 

mo = lim mefr(i) where mefr(t) = log{ ^^^ }^^ ) 

Note that since for the excited states rhi > m-o, then mefr(i) > rhesit + 1) > mo- This implies 
that the effective mass, defined above, is an upper bound on the ground state mass. In practice, 
sophisticated methods are used to choose loops G such that the correlation C{t) is dominated 
by the ground state glueball (ie to ensure |co| >> \ci\ ). By using a several different loops, 
a variational method can be used to achieve this effectively. These techniques are needed to 
obtain accurate estimates of mo from modest values of t since the signal to noise decreases as t 
is increased. Even so, it is worth keeping in mind that upper limits on the ground state mass 
are obtained in principle. 

The method also needs to be tuned to take account of the many glueballs: with different 
J^^ and different momenta. On the lattice the Lorentz symmetry is reduced to that of a 
hypercube. Non-zero momentum sates can be created (momentum is discrete in units of 2tt/L 
where L is the lattice spatial size). The usual relationship between energy and momentum is 
found for sufficiently small lattice spacing. Here we shall concentrate on the simplest case of 
zero momentum (obtained by summing the correlations over the whole spatial volume). 

For a state at rest, the rotational symmetry becomes a cubic symmetry. The lattice states 
(the Qi above) will transform under irreducible representations of this cubic symmetry group 
(called Oh)- These irreducible representations can be linked to the representations of the full 
rotation group SU(2). Thus, for example, the five spin components of a J^*-^ = 2"'"+ state 
should be appear as the two-dimensional E"^"*" and the three-dimensional Tg^^ representations 
on the lattice, with degenerate masses. This degeneracy requirement then provides a test for 
the restoration of rotational invariance - which is expected to occur at sufficiently small lattice 
spacing. 

The results of lattice measurements [||, |2|, |3|, of the 0^"'" and 2+"'" states are shown in fig 1. 
The restoration of rotational invariance is shown by the degeneracy of the two Oh representations 
that make up the 2~^^ state. Fig 1 shows the dimensionless combination of the lattice glueball 
mass m-o to a lattice quantity fo. We will return to describe the lattice determination of tq in 
more detail - here it suffices to accept it as a well measured quantity on the lattice that can be 
used to calibrate the lattice spacing and so explore the continuum limit. The quantity plotted, 
m-o^O) is expected to be equal to the product of continuum quantities moro up to corrections of 
order a^. This is indeed seen to be the case. The extrapolation to the continuum limit (a — > 0) 
can then be made with confidence. 

The value of ro in physical units is about 0.5 fm and we will adopt a scale equivalent to 
r^^ = 0.372 GeV. This information yields lattice predictions for the glueball masses of around 
1.6 GeV and 2.2 GeV for the 0"'"^ and 2"'""'" glueballs respectively. 

We return briefiy to the independence of the results on the volume of the lattice. In the 
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early days of glueball mass determination, it was expected that a spatial size L should satisfy 
m(0"'""'")L > 1 and, hence, that values of m(0"^^)L of 1 to 4 would suffice. A careful lattice 
study [ p!^ ] showed that m{Q'^~^)L > 8 was required to obtain rotational invariance and a result 
independent of L. The results collected in fig 1 all satisfy this latter inequality so can be regarded 
as the infinite volume determination. 

The predictions for the other J^^ states are that they lie higher in mass and the present 
state of knowledge is summarised in fig 2. Remember that the lattice results are strictly upper 
limits. For the J^*^ values not shown, these upper limits are too weak to be of use. 

Since quark - antiquark mesons can only have certain J^^ values, it is of special interest to 
look for glueball with J^'^ values not allowed for such mesons: ,0"' , 1 ^,2^ , etc. Such 
spin-exotic states, often called "oddballs" , would not mix directly with quark - antiquark mesons. 
This would make them a very clear experimental signal of the underlying glue dynamics. Various 
glueball models (bag models, flux tube models, QCD sum-rule inspired models,..) gave different 
predictions for the presence of such oddballs (eg. 1 ^) at relatively low masses. The lattice mass 
spectra clarify these uncertainties but, unfortunately for experimentalists, do not indicate any 
low-lying oddball candidates. The lightest candidate is from the T2 spin combination. Such 
a state could correspond to an 2"' oddball. Another interpretation is also possible, however, 
namely that a non-exotic 3"' state is responsible (this choice of interpretation can be resolved in 
principle by finding the degenerate 5 or 7 states of a J = 2 or 3 meson). The overall conclusion 
at present is that there is no evidence for any oddballs of mass less than 3 GeV. 

Glueballs are defined in the quenched approximation and, hence, they do not decay into 
mesons since that would require quark - antiquark creation. It is, nevertheless, still possible to 
estimate the strength of the matrix element between a glueball and a pair of mesons within the 
quenched approximation. For the glueball to be a relatively narrow state, this matrix element 
must be small. A very preliminary attempt has been made to estimate the size of the coupling of 
the 0'*'^ glueball to two pseudoscalar mesons . A relatively small value is found. Further work 
needs to be done to investigate this in more detail, in particular to study the mixing between 
the glueball and 0^^ mesons since this mixing may be an important factor in the decay process. 

Another lattice study will become feasible soon. This is to study the glueball spectrum in 
full QCD vacua with sea quarks of mass m£). For large rriD, the result is just the quenched result 
described above. For rriD equal to the experimental light quark masses, the results should just 
reproduce the experimental meson spectrum - with the resultant uncertainty between glueball 
interpretations and other interpretations. The lattice enables these uncertainties to be resolved 
in principle: one obtains the spectrum for a range of values of m/j between these limiting cases, 
so mapping glueball states at large m/5 to the experimental spectrum at light m£). 
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0.00 0,05 0,10 

Figure 1: The value of mass of the J^'~^ = 0^^ and 2^+ glueball states from refs[Q, ^, |3|, ^ in 
units of tq. The T2 and E representations are shown by octagons and diamonds respectively. 
The straight line shows a fit describing the approach to the continuum limit as a — > 0. 
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Figure 2: The mass of the glueball states with quantum numbers J^*-" from ref|Q. The scale 
is set by \/a « 0.44 GeV which yields the right hand scale in GeV. The solid points represent 
mass determinations whereas the open points are upper limits. 
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3 Potentials between quarks 



A very straightforward quantity to determine from lattice simulation is the interquark potential 
in the limit of very heavy quarks (static limit). This potential is of direct physical interest 
because solving the Schrodinger equation in such a potential provides a good approximation to 
the T spectrum. It is also relevant to exploring both confinement and asymptotic freedom on a 
lattice. 

The basic route to the static potential is to evaluate the average W{R, T) in the vacuum 
samples of a rectangular closed loop of colour flux (a Wilson loop of size RxT). here the lattice 
quantities R and T are related to the physical distances R and T by R = R/a, etc where a is 
the lattice spacing which is not known explicitly. Then it can be shown that the required static 
potential in lattice units is given by 



r-.oo W{R,T) 

The limit of large T is need to separate the required potential from excited potentials. This limit 
can be made tractable in practice by using more complicated loops than the simple rectangular 
loop described above. 

A summary of results |p for the potential at large R is shown in fig 3. The result that the 
force dV/dR tends to a constant at large R (and thus V{R) continues to rise as R increases) is 
a manifestation of the confinement of heavy quarks (in the quenched approximation) . The force 
appears to approach a constant at large R. A simple parametrisation is traditional in this field: 

V{R) = Vo-^+KR 

where K (sometimes written a) is the string tension. The term e/R is referred to as the 
Coulombic part in analogy to the electromagnetic case. The equivalent relationship in terms of 
quantities defined on a lattice is 

V{R) = %- ^ +KR 
R 

where the lattice string tension K = Ka?. 

Since the string tension is given by the slope of V{R) against R as R ^ oo, this implies that 
some error will arise in determining K coming from the extrapolation of lattice data at finite 
R. A practical resolution is to define a value of R where the potential takes a certain form. The 
convention is to use rg where 



^,dV{R) 



dR 



= 1.65 
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Figure 3: The force between static quarks at separation r as measured on lattice with different 
lattice spacing a (with the /3-values shown). 



Thus fo can be determined by interpolation in R rather than extrapolation. In practice, this 
means that fo is very accurately determined by lattice measurements and so is a useful quantity 
to use to set the scale since fg = ro/a. With the simple parametrisation above, we have 
Vq = (1.65 — e)/K so rg is closely related to the string tension since e ~ 0.25. The string 
tension is usually taken from experiment as ^/K = 0.44GeV where the value comes from cc and 
bb spectroscopy and from the light meson spectrum interpreted as excitations of a relativistic 
string. Similar analyses also imply that ro ~ 0.5 fm. Here we use Vq = 0.372 GeV to be specific. 
Since we shall be describing quenched lattice results, the energy scale set from different physical 
quantities will not necessarily agree (since experiment has full QCD not the quenched vacuum) 
and so a systematic error must be applied to any such choice of scale. This, for instance, must 
be kept in mind when taking glueball mass values from the lattice. 

The lattice potential V{R) can be used to determine the spectrum of bb mesons by solving 
Schrodinger's equation since the motion is reasonably approximated as non-relativistic. The lat- 
tice result is similar to the experimental T spectrum. The main difference is that the Coulombic 



part (e) is effectively too small (0.25 rather than 0.4). This produces |14] a ratio of mass dif- 
ferences (IP — 15')/(2S' — IS) of 0.71 to be compared with the experimental ratio of 0.78. This 
difference is understandable as a consequence of the Coulombic force at short distances which 
would be increased by 33/ (33 — 2Nf) in perturbation theory in full QCD compared to quenched 
QCD. We will return to discuss this. 
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Another feature of the lattice determination of potentials is that the energy of static quarks 
at separation R with an excited gluonic field can be determined. This enables predictions to be 
made for hybrid mesons (eg. bGb where G stands for the gluonic excitation). Such mesons can 
have J^^ values not allowed to bb states and this allows them to be explored experimentally. 
The current situation [|l5|, |l^ is that such states, as predicted by the lattice, will be at high 
masses and hard to isolate experimentally. 

At small R, the static potential can be used, in principle, to study the running coupling con- 
stant. Small R corresponds to large momentum and thus the coupling should decrease at small 
R. Thus the Coulombic coefficient e introduced above should actually decrease logarithmically 
as R decreases. Perturbation theory can be used to determine this behaviour of the potential 
at small R. 

In the continuum the potential between static quarks is known perturbatively to two loops 
in terms of the scale Aj^. For SU{3) colour, the continuum force is given by [Q] 

dV _ A a{R) 

dR ~ 3~R^ ^ ^ 

with the effective coupling a{R) given by 



47r[6olog(i?A^)-2 + (61/60) log log(i?AR)-2] 



(2) 



where 60 = ll/lGvr^ and 61 = 102 6o/121 are the usual coefficients in the perturbative expression 
for the /3-function, neglecting quark loops in the vacuum. Here Aji = 1.048Aj^. 

On a lattice the force can be estimated by a finite difference and one can extract the running 
coupling constant by using |Q 

* 2 ' 4 ' - Ri - Rj ^ ' 

where the error in using a finite difference is negligible in practice. 



This is plotted in fig 4 versus RVK: this combination is dimensionless and so can be deter- 
mined from lattice results since R^/K = Rv^, where K is taken from the fit to V{R)- The 
interpretation of a as defined above as an effective running coupling constant is only justified at 
small R where the perturbative expression dominates. Also shown are the two-loop perturbative 
results for a{R) for different values of A^. 

The figure clearly shows a running coupling constant. Moreover the result is consistent with 
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Figure 4: The effective running coupling constant a{R) obtained from the force between static 
quarks at separation R from refQ. The scale is set by the string tension K. Data are at /3 = 6.5 
(diamonds) and at /? = 6.2 (triangles). The dotted error bars represent an estimate of the 
systematic error due to lattice artefact correction. The curves are the two-loop perturbative 
expression with a(6.5)Aij = 0.060 (dotted) and 0.070 (continuous). 
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the expected perturbative dependence on R at small R. There are systematic errors, however. At 
larger R, the perturbative two-loop expression will not be an accurate estimate of the measured 
potentials, while at smaller R, the lattice artefact corrections (which arise because R/a ~ 1) 
are relatively big. Setting the scale using s/K = 0.44 GeV implies l/a(/? = 6.5) = 4.13 GeV, so 
R < 4a(6.5) corresponds to values of 1/i? > 1 GeV. This i?-region is expected to be adequately 
described by perturbation theory. 

This determination from the interquark force of the coupling a allows us to compare the result 
with the bare lattice coupling determined from /3 = G/g"^. At /3 = 6.5, Obarc = 5'^/4'7r = 0.073. 
the values of a shown in fig 4 are much larger. The effective coupling constant is thus almost 
twice the bare coupling. This is quite acceptable in a renormalisable field theory. The message 
is that the bare coupling should be disregarded - it is not a good expansion parameter. The 
measured a, however, proves to be a reasonable expansion parameter in the sense that the first 
few terms of the perturbation series converge. This successful calibration of perturbation theory 
on a lattice is important in practice. For instance, when matrix elements are measured on a 
lattice they have finite correction factors (usually called Z) to relate them to continuum matrix 
elements. These Z factors arc evaluated perturbatively - so an accurate continuum prediction 
needs trustworthy perturbative calculations. 

4 Full QCD 

So far we have discussed the glueball spectrum, interquark potentials and as in the quenched 
approximation. This corresponds to treating the sea quarks as of infinite mass (so they don't 
contribute to the vacuum). To make direct comparison with experiment, it is necessary to 
estimate the corrections from these dynamical quark loops in the vacuum. 

The strategy is to use a finite sea-quark mass but still a value larger than the empirical 
light quark mass. The reason is computational: the algorithms become very inefficient as the 
sea-quark mass is reduced. The target is to study the effects as the sea-quark mass is reduced 
and then extrapolate to the physical value. The present situation, in broad terms, is that there 
is no significant change as the sea-quark mass is reduced. This could be because there are no 
corrections to the quenched approximation. Alternatively, the corrections may only turn on at 
a much lower quark mass than has been explored so far. 

Let us try to make this argument a little more quantitative. For heavy sea-quarks of mass m, 
their contribution will be approximately proportional to e~^™/^ where E is a. typical hadronic 
energy scale (a few hundred MeV). Thus the quark loop contributions will be negligible for 
m >> E, which corresponds to the quenched approximation. As m ^ E, the effects will turn 
on in a non-linear way. 

The computational overhead of full QCD on a lattice is so large because the quark loops effec- 
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tively introduce a long range interaction. The quark interaction in the Lagrangian is quadratic 
and so can be integrated out analytically. This leaves an effective Lagrangian for the gluonic 
fields which couples together the fields at all sites. This implies that, in a Monte Carlo method, 
a change in gluon field at one site involves the evaluation of the interaction with all other sites. 
In practice, one makes small changes at all sites in parallel, but this still amounts to inverting 
a large sparse matrix for each update. This is computationally slow. 

As the sea-quark mass becomes small, one would expect to need a larger lattice size to hold 
the quarks. For heavy quarks, the effective range of the quark loops in the vacuum will be of 
order 1/m. Thus the quenched approximation corresponds to m — > oo and a local interaction. 
For light quarks of a few MeV mass, the range will not be 1/m, because quarks are confined. The 
lattice studies that have been made suggest that spatial sizes of order twice those adequate for the 
quenched approximation are needed for full QCD. This also implies considerable computational 
commitment. 

A popular indication of how close a full QCD study is to experiment is to ask whether the 
p meson can decay to two pions. Since the decay is P-wave, it needs non-zero momentum. On 
a lattice spatial momentum is quantised in units of 271 /L. Thus we need my > 2mp -\- 2tt/L 
for the decay channel to be allowed energetically. At present this criterion is rarely satisfied in 
quenched studies, let alone in the more computationally demanding case of full QCD. 

The conclusion of current full QCD lattice calculations is that the expected sea-quark effects 
are not yet fully present. The main effect observed in full QCD calculations is that the lattice 
parameter (3 which multiples the gluonic interaction term in the Lagrangian is shifted. Apart 
from this renormalisation of /?, there is little sign of any other statistically significant non- 
perturbative effect. 

Consider the changes to be expected for the inter-quark potential when the full QCD vacuum 
is used: 



• At small separation i2, the quark loops will increase the size of the effective coupling a 
compared to the pure gluonic case. This effect can be estimated in perturbation theory 
and the change at lowest order will be from 1/33 to l/(33-2A'^^). 

• At large separation R, the potential energy will saturate at a value corresponding to two 
'heavy-quark mesons'. In other words, the flux tube between the static quarks will break 
by the creation of a qq pair from the vacuum. 



Current lattice simulation [13| shows some evidence for the former effect but no statistically 



significant signal for the latter. 

If one assumes that these lattice simulations are an approximation to the true full QCD 
vacuum, then one can use them to estimate the full QCD running coupling from lattice studies. 
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A summary [|To| of present lattice results is that 

«M-s(^^) = 0-112(7) 

This conclusion will be reinforced when the full QCD lattice results reproduce more of the 
features of the experimental spectrum. 



5 Quenched hadron masses 

Since the full QCD simulation is not feasible at present, it is worthwhile to explore the hadron 
spectrum in the quenched approximation. This amounts to allowing quarks to propagate in the 
gluonic vacuum. Computationally this can be studied by solving the lattice Dirac equation for 
the quarks. Since the gluonic vacuum is full of rich structure, this is a computationally intensive 
problem: it amounts to inverting a large sparse matrix. Indeed it proves necessary to compute 
the quark propagation for a range of valence quark masses larger than the physical light quark 
masses and then extrapolate. Because of this, the statistical precision of such calculations is 
still somewhat limited. Nevertheless, different groups using different methods agree on the main 



results. For a recent review see ref[10|. 



Let us first discuss the p meson. In fig 5, the mass is compared [16| to a lattice scale (here 
E = Ka? where K is the string tension). The figure shows that different treatments of fermions 
on the lattice (Wilson, Clover and staggered) with different discretisation errors (a, a^o and 
a? respectively) are in agreement with a common continuum limit. Moreover, using the usual 
convention that the string tension is ^/~K = 0.44 GeV, the continuum value of mp/ \fK of 1.80(5) 
yields rrip = 0.79(2) GeV which is consistent with the experimental value of 0.77 GeV. 

Studies have been made of the mesons and baryons which are composed of light and strange 



quarks- see ref[10|. The surprise is that the quenched approximation seems to reproduce these 
experimental values quite well. This may reflect the relatively large errors that are still present 
in the lattice determinations. It may also reflect the fact that the hadronic dynamics has a 
similar energy scale in each case so that the quenched approximation makes similar errors - 
which cancel in mass ratios. 



6 Matrix elements eg fp 

One of the advantages of lattice QCD is that it is a method to calculate hadronic matrix elements 
from first principles. Consider as an example the weak decay of a pseudoscalar meson P. The 
weak axial current will couple to the quarks in the meson. This current will be local. Thus the 
required quantity will relate the quark current to the hadronic state P. This is the pseudoscalar 
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Figure 5: The mass of the p meson as a dimensionless ratio to the lattice string tension E = Ka^. 
The continuum limit corresponds to a ^ at the left hand side. The filled squares are with 
staggered fermions, the open squares from Wilson fermions and the triangles are from Clover 
fermions. The lines show the discretisation errors which behave as a^, a and a^a respectively. 
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Figure 6: The diagram representing a static quark Q propagating in the time direction and a 
hght quark q which creates a gauge invariant quantity appropriate for studying /b- 

decay constant fp which is defined by the matrix element of the divergence of the axial current. 
For a pseudoscalar state of zero spatial momentum, 



where j is the axial current which in terms of quark fields is ^757^9. For the pion, this identity 
is the partially conserved axial current (PCAC) relation and /^r is the coupling of the pion to 
the weak current (and hence is relevant to the fx + v decay mode of the pion) . Since the axial 
current is represented by local quark fields, fp gives a relationship between the hadronic state 
{\P >) and the quark sub-structure (j^)- 

For a pseudoscalar meson with a heavy quark, such as the B meson, this same relationship is 
needed. Because the relevant weak decay is not currently observable (branching ratio too small), 
a lattice calculation is needed to determine /b- The picture, for a heavy quark, is clarified by 
the heavy quark effective theory (HQET) which treats the heavy quark as slow moving. In 
the static limit, the lattice calculation needed is easily visualised. A straight line of colour flux 
of length T in the lattice time direction represents the propagation of the heavy quark. By 
combining this with the propagator from one end to the other of a light quark in the lattice 
gluonic background field, one has a gauge invariant quantity which can be measured on a lattice 
(see fig 6). At each end, one joins the heavy quark to the light quark by a local axial current 
(7574 for a state at rest). Then the observed correlation C{t) is proportional to fpe~^P^ for 



fp 



\J2mp 
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large T so allowing fp and Ep do be determined in principle. In practice more sophisticated 
methods are used to improve the lattice measurement signal. Here Ep has the interpretation 
of the mass difference of the B meson and the h quark - although this difference is not directly 
useful since a non-perturbative definition of the quark mass requires careful discussion. 

As well as using static quarks, lattice studies have been performed using propagating heavy 
quarks. Combining the results from both studies ||l^ yields = 180(50) MeV. This value is 
needed as an ingredient in using experimental data to fix the CKM weak matrix elements. This 
is turn has implications for the experimental feasibility of CP violation through studies of BB 
mixing. 

7 Outlook 

Lattice techniques can extract reliable continuum properties from QCD. At present, the com- 
putational power available combined with the best algorithms suffices to give accurate results 
for many quantities in the quenched approximation. The future is to establish accurate values 
for more subtle quantities in the quenched approximation (eg. weak matrix elements of strange 
particles) and to establish the validity of the quenched approximation by full QCD calculations. 

We need to reach the stage where an experimentalist saying 'as calculated in QCD' is assumed 
to be speaking of non-perturbative lattice calculations rather than perturbative estimates only. 
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Scalar glueball decay* 



J. Sexton^, A. Vaccarino and D. Weingarten 
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We evaluate the coupling constant for the lightest scalar glueball to decay to pseudoscalar meson pairs. The 
calculation is done in the valence approximation on a 16'^ x 24 lattice at /? = 5.70 for two different values of 
pseudoscalar meson mass. 



A valence approximation calculation [1] using 
between 25000 and 30000 gauge configurations 
at each of several different [3 gives a value of 
1740 ±71 MeV for the infinite volume contin- 
uum limit of the mass of the lightest scalar glue- 
ball. An independent calculation [2] with between 
1000 and 3000 configurations at several different 
[3, when interpreted [3] following Ref. [1], yields 
1625 ± 94 MeV for the lightest scalar glueball 
mass. Finally, an analysis in Ref. [3] suggests 
that the error in the valence approximation to 
the scalar glueball mass is probably less than 60 
MeV. The mass calculation with larger statistics 
then appears to favor /j(1710) as a scalar glue- 
ball, with /o(1590) a possible but less likely can- 
didate. A third alternative is that the scalar glue- 
ball is a linear combination of the two states with 
significant contributions from each. No other well 
established resonance consistent with scalar glue- 
ball quantum numbers lies within 200 MeV of the 
1740 MeV prediction. 

The interpretation of the scalar glueball as 
/j(1710) or some combination of /j(1710) and 
/o(1590) is obscured, however, by two unresolved 
issues. First, if the scalar glueball decay width 
were large enough, this particle would not be eas- 
ily identified, in which case it need not be any well 
established resonance. Second, since the scalar 
glueball is a fiavor singlet, one might expect its 
coupling constant to a pair of rj' s to be the same 
as its coupling to a pair of K' s or a pair of tt's. In 
fact, the observed values of these three coupling 
constants are far from equal both in the case of 
/j(1710) and in the case of /o(1590). 

We are now in the process of calculating 
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the scalar glueball's coupling to pairs of pseu- 
doscalars. Our preliminary results strongly sug- 
gest this particle has a width significantly below 
200 MeV. Thus we consider it quite likely that 
the scalar glueball has already been found in ex- 
periment. In addition, the coupling constants we 
obtain show a variation with the mass of the de- 
cay products consistent with the variation in both 
the /j(1710) and the /o(1590) couplings. 

To evaluate the glueball decay couplings, con- 
sider a euclidean lattice gauge theory, on a lattice 
X T, with the plaquette action for the gauge 
field, the Wilson action for quarks, and the gauge 
field transformed to lattice Coulomb gauge. Let 
U[(x) for a space direction i = 1,2,3, and integer 
radius r > 1, be a smeared link field found by 
averaging the links running in direction i from 
the sites of the (r — 1) x (r — 1) square oriented 
in the two positive space directions orthogonal to 
i starting at site x. Let Vfj(x) be the trace of the 
product of U[(y) and Uj (y) around the boundary 
of an (r — 1) X (r— 1) square in the ij plane begin- 
ning at X. Let the zero-momentum scalar glueball 
operator g''(t) be the sum of the Vfj(x) for all i,j 
and X in the time t lattice hyperplane. Define the 

r' I 

quark and antiquark fields ^ {x) and {x) to 
be Wilson quark and antiquark fields smeared by 
convoluting the local Wilson fields with a space 
direction gaussian with root-mean-square radius 
r'. Define the smeared pseudoscalar field 7r[ {x) 

with fiavor index i to be ^ (a;)7®Aj'^'' {x), where 
\i is a Gell-Mann fiavor matrix. Let 7r[ {k, t) be 
the Fourier transform of 7r[ {x) on the time t lat- 
tice hyperplane. We now fix r and r' to the values 
3 and \/&, respectively, used in the summary of 
results presented below, and omit the radius su- 
perscripts. 
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At large time separations t, the vacuum ex- 
pectation values < ^1(0, t)^j(0, 0) > and < 
TTj(k,t)TTi(k,0) > approach ri'fL^exp(—E,^it) and 

'ri2L^exp(—ET^2i), respectively. Here k is any 
momentum with norm 2'itL~^, Ej^i is the mass 
of a pion, and £"^2 is the energy of a pion 
with momentum 2'itL~^. In addition, in the va- 
lence approximation, < g(t)g(0)) > approaches 
'ri'^L^exp(—Egt), where Eg is the mass of the light- 
est scalar glueball. At this point it is convenient 
to redefine g(t), TTi(0,t) and TTi(k,t) by dividing 
each by the corresponding rj. 

Define the two pion zero momentum opera- 
tor tl(0,t) to be 6-1/2^. *i(0,t)*i(0,0). Let 
the two pion operator Tl(k,t) for nonzero k be 
6~^^.^Tti(k,t)Tri(—k,0) where the sum is over 
the six ways of placing a vector k of length k 
in a positive or negative lattice coordinate direc- 
tion. Define the vacumm subtracted operators 
f[s(k,t) to be fl(k,t)- < n\fl(k,t)\n >. Let |1 > 
be the lowest energy eigenstate in ns(0,t)|Sl >, 
and let |2 > be the lowest energy eigenstate in 
fls(27r_L-i,t)|fi >. The states |1 > and |2 > are 
both normalized to 1. Let £',r7ri and £',r7r2 be the 
energies of these two states, respectively. Define 
the amplitudes riij(t) 

r,ii{t) = 1.-3 < i|fl,(0,O|fi>, (1) 
rii2{i) = L-"" <i\Il,{27rL-\t)\n>, (2) 

For large t, riij(t) has the asymptotic form 
riijexp(—ET^jt). Connected three-point functions 
from which coupling constants can be extracted 
are now given by 



<g(tg)fis(0,U)>, (3) 
< g{tg)Il,{27rL-\t,) > . (4) 



If the quark mass, and thus the pion mass, is 
chosen so that £',r7ri is equal to Eg, the lightest 
intermediate state which can appear between the 
glueball and pions in a transfer matrix expression 
for Ti(tg,tT^) is |1 >. Thus for large enough tg 
with fixed, Ti(tg,t^) will be proportional to 
the coupling contant of a glueball to two pions at 
rest. If the quark mass is chosen so that £',r7r2 
is equal to Eg, however, the lightest intermediate 
state which can appear between the glueball and 



pions in a transfer matrix expression for 72 (t^ , t,r) 
is still |1 >, not |2 >. This holds since ?yi2(0 
in Eq. (2) has no reason to vanish. To obtain 
from T2(tg,tT^) the coupling of a glueball to two 
pions with momentum 2'itL~^, the contribution 
to T2(tg,tT^) arising from the |1 > intermediate 
state must be cancelled off. From the three-point 
functions defined in Eqs. (3) and (4) we therefore 
define the amplitudes 

S,itg,t,) = mg,U) - (5) 

Tl22\l-K ) 

S2{tg,t,) = T2{tg,t,) - ^^^Ti(t,,tO. (6) 
Tliiyi-K ) 

In S2{tg,tTi) the contribution of the undesirable 
|1 > intermediate state has been cancelled. In 
Siiigjiw) a contribution from the intermediate 
state |2 > has been cancelled. Although the sub- 
traction in Si(tg, tir) is irrelevant for large enough 
tg, we expect that as a result of this subtraction 
Siiigjiw) will approach its large tg behavior more 
rapidly than does Ti(tg,t^). 

An additional intermediate state which can also 
appear in a transfer matrix expression for either 
Ti(tg,tT^) is the isosinglet scalar bound state of a 
quark and an antiquark. The relative contribu- 
tion of this state to either three-point function, 
however, is suppressed by the factor iJ,\tg — 
where /j, is the coupling constant of this state to 
the glueball. Phenomenological arguments fol- 
lowing Ref. [3] strongly suggest that /j, is less than 
60 MeV. In our main results below, removing the 
additional intermediate state would therefore give 
at most a 4% correction. 

At large tg and , the three-point functions 
become 



Sl{tg,t,r) 

V3Airyii(l - ri2r2i)L^ 



S2{'tg , iw) 



3A2T?22(1 - ri2r2i)L^ 

V^e;ei2 



S2itg,tn)- 



(7) 



(8) 



Here Ai and A2 are the glueball coupling con- 
stants to pseudoscalar pairs at rest or with mo- 
mentum 2'itL~^, respectively. Up to a factor of 



2 



—i these coupling constants are invariant decay 
amplitudes with the standard normalization con- 
vention. The factors rjij are given by the large t 
behavior of riij(t) as discussed earlier, and rij is 
Tjij/rjjj. The factors Si(tg,tT^) are 

Si{tg,t,,) = '^exp{-Eg\t -tg\- 

t 

E,i\t-U\- E,i\t\). (9) 

In Eq. (9) we make the simplifying assumption 
that -E'lTTri is 'iEiri. In our results to be presented 
below we actually use a slightly more complicated 
expression for Si(tg,tT^) which takes into account 
small differences between E'jrTri and 2£'^j'. 

To obtain values of A,- from Eqs. (7) and (8) 
we need the amplitudes riij(t). These we deter- 
mine from propagators for two pion states. Define 
the two pion operator Tl(ti , t'^) from smeared pion 
fields, not Fourier transformed and not renormal- 
ized, to be 6""'^''^ 7rj(0, ti)7rj(0, ^2)- Define two 
pion propagators by 

Ci{tiM) = 

<ii{h + 2h,h + h)nsio,t2)>, (10) 
^2(^1,^2) = 

< n{ti+2t2,ti+t2)tl,{27r/L,t2) > . (11) 
For large values of ti, these amplitudes approach 

Ci{ti , ^2) = 

Ciiexp(-ET,T,iti) + C2iexp(-ET,T,2ti), (12) 
dj = 

ViVij(i2)Vii(i2) + VQV2Vij(i2)Vi2(i2), (13) 

where the rji are determined from the asymptotic 
behavior of single pion propagators as discussed 
earlier. 

At l3 of 5.7, we found that k of 0.1650 nearly 
gives -E'ttttI equal to Eg, and k of 0.1675 nearly 
gives £',r7r2 equal to Eg. For k of 0.1650 and t from 
to 5, we determined riij(t) with Eqs. (12) and 
(13) applied to Ci{ti,t2) from an ensemble of 100 
independent gauge configurations on a 16^ x 40 
lattice. For k of 0.1675 and t of to 5, we deter- 
mined riij(t) from Ci{ti,t2) using 870 indepen- 
dent gauge configuratons. The riij(t) turn out to 
be nearly independent of t. In each case ryoo is 
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Figure 1. Decay couplings. 

statistically consistent with 1.0, rjn is between 
1.0 and 1.1 and ryi2 and — ?y2i are positive and 
less than 0.1. 

We then evaluated Ai for k of 0.1650 and A2 
for k of 0.1675 using Eqs. (7) and (8) applied 
to Si(tg,tT^) and S2{tg,tTi) found with an ensem- 
ble of 7200 independent gauge configurations on 
a lattice 16^ x 24. We obtained statistically sig- 
nificant results for tg—tj^ of 0, 1 and 2 and of 1, 
2, 3, 4, and 5. The A,- are nearly constant, within 
statistical errors, across these 15 points. Figure 1 
shows predicted coupling constants fitted to the 
data with tg — of 1 in comparision to observed 
decay couplings for decays of /j(1710) to pairs of 
rj's, K's and it's. The horizontal axis gives the 
squared mass of the decay products. Masses and 
decay constants are shown in units of the p mass. 
If couplings to rj' s, K' s and tt's are determined 
from our predicted couplings by linear interpola- 
tion in the squared mass of the decay products, we 
predict a total width for glueball decay to pseu- 
doscalar pairs of 52t^^ MeV, in comparison to 
84± 23 MeV for /j(1710). 

We would like to thank Steve Sharpe for finding 
an error in an earlier version of this work. 
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